Paper No. T.1.3 




SHOCK/VORTEX INTERACTION AND 
VORTEX-BREAKDOWN MODES 


O. A. Kandil and H. A. Kandil 

Old Dominion University, Norfolk, VA, USA 


C. H. Liu 

NASA Langley Research Center, Hampton, VA, USA 


»' Al (ac< , 






7 ®sa roK -' 0 


IUTAM SYMPOSIUM OF FLUID 
DYNAMICS OF HIGH ANGLE OF ATTACK 
Tokyo, Japan — September 13-17, 1992 





SHOCK/VORTEX INTERACTION AND 
VORTEX-BREAKDOWN MODES 


O. A. Kandil and H. A. Kandil 
Dept, of Mechanical Engineering and Mechanics 
Old Dominion University, Norfolk, VA 

, C. H. Liu 
Theoretical Flow Physics Br. 

NASA Langley, Hampton, VA 


Summary 

Computational simulation and study of shock/vortex interaction and vortex-breakdown modes 
are considered for bound (internal) and unbound (external) flow domains. The problem is 
formulated using the unsteady, compressible, full Navier-Stokes (NS) equations which are 
solved using an implicit, flux-difference splitting, finite-volume scheme. For the bound flow 
domain, a supersonic swirling flow is considered in a configured circular duct and the problem 
is solved for quasi-axisymmetric and three-dimensional flows. For the unbound domain, 
a supersonic swirling flow issued from a nozzle into a uniform supersonic flow of lower 
Mach number is considered for quasi-axisymmetric and three-dimensional flows. The results 
show several modes of breakdown; e.g., no-breakdown, transient single-bubble breakdown, 
transient multi-bubble breakdown, periodic multi-bubble multi-frequency breakdown and 
helical breakdown. 

Introduction 


Longitudinal vortex/transverse shock-wave interactions are typical applications which appear 
in transonic and supersonic flows over a strake-wing configuration at moderate-high angles 
of attack, at a supersonic inlet injesting a vortex and inside a supersonic combustor where 
fuel is injected in a swirling jet to enhance fuel-air mixing [l]-[3]. For the strake-wing 
configuration, vortex breakdown is undesirable since it results in the stall phenomenon, and 
hence its occurrence need to be delayed. On the other hand, vortex breakdown for the other 
two applications is desirable since it enhances mixing and stability of the flame [4]-[5], and 
hence its occurrence need to be controlled for the optimum performance. Unfortunately, the 
literature lacks this type of analysis with the exception of the preliminary work of Liu, Krause 
and Menne [6], Copening and Anderson [7], Delery, et al. [1], Kandil and Kandil [8] and 
Meadows, Kumer and Hussaini [9]. 

The first time-accurate NS solution for a supersonic vortex breakdown was developed by 
the present authors in Ref. [10]. We considered a supersonic quasi-axisymmetric vortex 
flow in a configured circular duct The time-accurate solution of the unsteady, compressible 
NS equations was obtained using an implicit, upwind, flux-difference splitting finite-volume 
scheme. A shock wave has been generated near the duct inlet and unsteady vortex breakdown 
has been predicted behind the shock. The predicted flow was characterized by the evolution, 
convection merging and shedding of vortex breakdown bubbles. The Euler equations were 



also used to solve the same problem. The Euler solution showed larger size and number of 
vortex-breakdown bubbles in comparison with those of the NS solutions. The time-accurate 
solution was carried out for 3,200 times steps which are equivalent to a dimensionless time 
of 16. Only one value of Reynolds number of 10,000 was considered in Ref. [10]. 

In a later paper [11], we expanded our study of this flow using time-accurate computations 
of the NS equations with a fine grid in the shock- vortex interaction region and for longer 
computational times. Several issues were addressed in that study. First, we showed the 
effect of Reynolds number on the temporal evolution and persistence of vortex- breakdown 
bubbles behind the shock. In that stage of computations, the conditions at the downstream 
exit were obtained by extrapolating the components of the flowfield vector from the interior 
cell centers. Although the flow was supersonic over a large portion of the duct exit, subsonic 
flow existed over a small portion of the exit around the duct centerline. Therefore, selected 
flow cases were computed using a Riemann-invariant-type boundary conditions as well as 
other boundary conditions at subsonic points of the duct exit [12]. 

In the present paper, we consider shock/vortex interaction and the resulting vortex breakdown 
modes for quasi-axisymmetric and three-dimensional flows. This study covers bound and 
unbound flow domains. For the bound domain, supersonic swirling flow is considered in a 
configured duct, and for the unbound domain supersonic swirling flow that is issued from a 
nozzle into a uniform supersonic flow of lower Mach number is considered. 

Highlights of Formulation and Computational Scheme 


Formulation: The conservative, unsteady, compressible, full Navier-Stokes equations in 

terms of time-independent, body-conformed coordinates £*, and are used to solve the 
problem. The equations are given in Ref. [11] and hence they are not presented here. Along 
with these equations, boundary conditions are specified at the computational-domain inlet, 
side wall and downstream exit. The downstream exit boundary conditions will be presented 
and discussed in the next section of the computational results. The initial conditions are also 
presented in the next section. 

Computational Scheme: The computational scheme used to solve the unsteady, compress- 
ible full NS equations is an implicit, upwind, flux-difference splitting, finite-volume scheme. 
It employs the flux-difference splitting scheme of Roe which is based on the solution of the 
approximate one-dimensional Riemann problem in each of the three directions. In the Roe 
scheme, the inviscid flux difference at the interface of a computational cell is split into left 
and right flux differences. The splitting is accomplished according to the signs of the eigen- 
values of the Roe averaged-Jacobian matrix of the inviscid flux at the cell interface. The 
smooth limiter is used to eliminate oscillations in the shock region. The viscous and heat-flux 
terms are differenced using second-order spatially accurate central differencing. The resulting 


difference equation is approximately factored and is solved in three sweeps in the £ 2 and 
£ 3 directions. The scheme is coded in the computer program which is called “FTNS3D”. 

The quasi-axisymmetric solutions are obtained using the three-dimensional code by forcing 
the flowfield vector to be equal on two axial planes, which are in close proximity of each other. 

Computational Results and Discussion 


a, Quasi-axisymmetric Vortex Breakdown Modes in a Configured Duct: Figure 1 shows 
an axial plane of a configured circular duct. The design of the duct configuration is intended to 
ensure that the supersonic inflow becomes supersonic at the exit. However, as the computation 
will show, a small portion of the duct exit becomes subsonic at certain times for the specified 
inflow conditions. A grid of 221x51x2 in the axial, radial and circumferential directions, 
respectively, is used for the quasi-axisymmetric flow. For the three-dimensional flow, a grid 
of 221x51x49 in the axial, radial and circumferential directions, respectively, is used. The 
minimum grid length is 0.002. The upstream Mach number is kept at 1.75. 


The initial profile for the tangential velocity is given by 


w 




1 — exp 



( 1 ) 


where Uoo = 1.74, r m = 0.2 and k e = 0.1. The maximum p^, swirl ratio 0, is at r = 
0.224 and its value is kept at 0.32. The radial velocity, v, at the initial station is set equal 
to zero and the radial momentum equation is integrated to obtain the initial pressure profile. 
Finally, the density p is obtained from the definition of the speed of sound for the inlet flow. 
With these compatible set of profiles, the computations are carried out accurately in time with 
A t - 0.0025. The wall boundary conditions follow the typical Navier-Stokes solid-boundary 
conditions. 


Reynolds number, Re = 4,000: Figure 2 shows snapshots of the streamlines and Mach 

contours for the flow case of Re = 4,000. For this value of Reynolds number, a single 
breakdown bubble is seen at t = 5 and it is convected downstreams as time passes. This 
breakdown bubble is formed during the downstream motion of the inlet shock, which reaches 
its maximum downstream displacement at t = 5. Later on, the shock moves upstream, as 
it is seen at t = 8, while the breakdown bubble is convected in the downstream direction. 
Thereafter, the shock stays stationary at the inlet This swirling flow case shows a transient 
single-bubble breakdown flow. The conditions at the exit are obtained by extrapolation from 
the interior cell centers. 

Reynolds number, Re = 20,000: Figure 3 shows snapshots of the streamlines and Mach 

contours for this case. These snapshots show a vortex breakdown mechanism of evolution. 




convection, merging and shedding of bubbles while the inlet shock is moving downstreams, 
then upstreams and finally downstreams. The inlet shock becomes stationary and no more 
bubbles are developed This swirling flow case shows a transient multi-bubble breakdown 
flow. 

Reynolds number, Re = 100,000: Figure 4 shows snapshots of the streamlines and Mach 

contours for this case. The downstream boundary conditions are obtained by extrapolating 
all the flow variables from the cell centers at the exit. The streamlines snapshots show multi- 
bubble vortex breakdown evolution, convection, merging and shedding. The time-accurate 
integration was carried out up to t = 200 and the solution showed periodic multi-frequency 
cycles of vortex-breakdown bubbles [21]. An example of the merging of vortex breakdown 
bubbles of same sign of vorticity is shown at t = 17. An example of convection and shedding 
of vortex breakdown bubbles is shown at t = 25. Comparing the streamlines solutions at t 
= 25 and t = 89, it is seen that the solutions are almost the same which conclusively show 
that the breakdown process is periodic. The Mach-contours show the dynamics of inlet shock 
motion. In the time range of t = 3-8, the inlet shock moves upstream toward the inlet and 
its central portion exists outside the inlet section at t = 8. In the time range of t = 8-25, 
the inlet shock moves downstream with corresponding evolution, convection, merging and 
shedding of breakdown bubbles. In the time range of t = 25-45, the inlet shock maintains its 
motion in the downstream direction at a slower rate than before, while another shock, which 
is downstream of the inlet shock, appears and also moves in the downstream direction. The 
evolution, convection and shedding slowly continues until t = 66. In the time range of t = 66- 
78, the downstream shock disappears and a large vortex-breakdown bubble appears and moves 
upstream. This motion of the bubble is accompanied by upstream motion of the inlet shock 
(t = 78). Later the inlet shock again moves in the downstream direction and the process is 
repeated. An animation movie has been produced for the total dimensionless time of t = 200. 

Exit Riemann Invariant Condition, = 2poo: In this case, the back pressure at the 

subsonic points of the duct exit is specified to be p* = 2poo and the other four variables 
are extrapolated from the interior cell centers. Figure 5 shows snapshots of the streamlines 
and Mach contours of the solution. Comparing the present solution with the solution of the 
previous case (Fig. 4), it is seen that the two solutions are similar with the exception that the 
present solution lags that of the first case in time. The reason behind this behavior is that the 
back pressure pj is larger than that of the first case. Moreover, the Riemann invariant type 
conditions at subsonic points let the downstream effects propagate upstream as time increases. 
The existence of the large back pressure which is felt upstream supports the inlet shock and 
keeps it in the . inlet region. 

b. Three-Dimensional Vortex-Breakdown in a Configured Duct: Figure 6 shows the 

streamlines and Mach contours for the three-dimensional vortex breakdown modes in a con- 


figured duct with the same inflow conditions as those of Fig. 4. At the downstream exit, flow 
conditions are extrapolated from the interior cell centers. It is observed from the vertical, 
axial planes that the vortex breakdown is a three-dimensional phenomenon. The breakdown 
changes from a two-bubble type (t = 2-5) to a multi-bubble type (t = 7-34). At certain 
times, t = 11 and 19, it is observed that the multi-bubble breakdown is followed by a spiral 
breakdown. Beyond t = 46, the breakdown becomes a single-bubble type. 

c. Quasi-axisymmetric Vortex-Breakdown of a Supersonic Flow from a Nozzle: Here, a 
supersonic swirling jet at = 3.0, which is issued from a nozzle into a supersonic uniform 
flow of = 2.0, is considered. A grid of 221x51x2 in the axial, radial and tangential 
directions, respectively is considered. The computational domain in an axial plane has the 
dimensions of 7.0 x 3.5 in the axial and radial directions, respectively, where the nozzle exit 
radius r = 1. The freestream Reynolds number is 296,000. The inflow profiles are taken 
from the experimental data of Ref. [2] and they are used as quasi-axisymmetric profiles [12]. 
Figure 7 shows snapshots of the streamlines and Mach contours at selected time steps. The 
streamline figures show multi-bubble breakdown at the early levels of time. These bubbles 
develop due to the shock system existing at the nozzle exit. The shock system consists of 
two oblique shocks; one is weak and the other is strong. For t > 55, the shock system and 
the breakdown bubble move upstreams and remain there with low frequency oscillation. 

d. Three-Dimensional Vortex-Breakdown of a Supersonic Flow from a Nozzle: The flow 
of the previous case is considered for three-dimensional computations using a square-cross 
flow computational domain and a circular cross-flow computational domain. For the first case, 
one-half the square side is 3.5 units and for the second case, the radius is 3.5 units. Three 
types of grid are used for the first case (Gird types 1, 2 and 3) and one type of grid is used for 
the second case (Grid type 4), see Fig. 8. This study shows the dominant effect of the grid 
shape and distribution on the vortex breakdown solutions. While grid type 1 shows two large 
breakdown bubbles, Grids 2, 3 and 4 show substantially smaller breakdown bubbles. More- 
over, the conical shock system ahead of the breakdown bubbles shows substantial difference 
in size, location and shape. This study gives an alarming conclusion to the researchers who 
use substantially coarse grids with a rectangular cross-flow domain and a rectangular grid, 
and still claim capturing of physical vortex-breakdown flows. In the present study, the vortex 
breakdown bubbles of grid type 4 are of the same size order as those of grid type 3, although 
the shock system is somewhat different in shape. 

Concluding Remarks 


Computational simulation and study of shock/vortex interaction and vortex breakdown have 
been considered for internal and external supersonic swirling flow. The time-accurate 
computation for full Navier-Stokes equations is used to produce all the present cases. 



Numerous modes of vortex breakdown flows have been captured. The results show the effects 
of Reynolds number, downstream exit conditions, and grid shape, fineness and distribution. 
Three-dimensional vortex breakdown computations show somewhat different results from the 
quasi-axisymmetric counter parts since the three-dimensional flow provides a flow reliefing 
effect in the circumferential direction. 
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Figure 1. Configured duct grid 


Figure 2. Streamlines and Mach contours for a swirling flow with a transient single-bubble 
breakdown. Moo = 1.75, /? = 0.32, Re = 4,000. 


t = 15 


r = 17 


t = 30 


Figure 3. Streamlines for a swirling flow with transient multi-bubble breakdown, - 
0 = 0.32, Re = 20,000 

































t = 2 



Figure 6. Mach contours and streamlines for a swirling three-dimensional flow in a configured 
circular duct. Moo = 1-75, Re = 100,000, 0 = 0.32 
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ABSTRACT 

The near-vortex-wake flow of a large aspect-ratio rect- 
angular wing is accurately computed by using the thin- 
layer and full Navier-Stokes (NS) equations. The chord- 
wise section of the wing is a NACA-0012 airfoil and its 
tip is round. The computations have been carried out on 
a fine C-0 grid using an implicit, upwind, flux-difference 
splitting, finite-volume scheme. The thin-layer NS results 
have been obtained with and without flux limiters and 
the full NS results have been obtained without flux lim- 
iters. Flow transition from laminar to turbulent is mim- 
iced by tuming-on the Baldwin-Lomax algebraic model 
at an experimentally prescribed chord-station location of 
0.05. Comparison of the computed results with each other 
and with the experimental data shows that the full NS 
results give the best resolution of the near-vortex-wake 
flow. Next, the strength of the wing-tip vortex has been 
reduced substantially without reducing the lift coefficient 
by using flow-injection from a slot along a portion of the 
wing tip. The flow injection is directed in the wing plane 
at 45° with the wing-tip chord. 

INTRODUCTION 

Recently, there has been an alarming rate of increase 
of the volume of air civil transport at many airports. With 
this status of air traffic congestion, the hazardous effect 
of wake vortices emanating from large aircraft on light 
aircraft, which pass through these vortices during land- 
ing and take-off operations, is becoming a vital concern 
for operations safety. These vortices are characterized 
with high intensity and turbulence, and may produce high 
rolling moments on trailing aircraft which could exceed 
the available roll control of the trailing aircraft More- 
over, the trailing aircraft, under the influence of these 
vortex trails, could suffer a loss of altitude or climb rate 
in addition to structural damages. These vortices may 
persist up to several miles and for long periods of time 
before their decay, and therefore they play a major role 
in sequencing landing and take-off operations at busy air- 
ports. For example, a minimum safe-separation distance 
between aircraft is dictated depending on the vortex inten- 
sity, wind shear, atmospheric turbulence and temperature 
gradient; among others. 
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The literature shows a few experimental and compu- 
tational investigations that attempt to model and analyze 
wake vortex interaction, merging, decay and their haz- 
ardous effects on trailing aircraft Hallock and Eberie 1 
presented a review of the state of the art of aircraft wake 
vortices covering the research efforts in the United States 
until the mid-seventies. Experimental wind tunnel and 
airport measurements of the vortex wakes were conducted 
by Dee and Nicholas 2 Harvey and Perry 3 , El-Ram ly 4 . 
Wood and McWilliams 5 , Gardoz 6 , Cliffone and Lonzo 7 , 
Olwi and Ghazi*, Liu 9 and Liu et al 10 . 

Mathematical models and computational schemes 
were developed using inviscid analysis by Chorin and 
Bernard 11 , Hacket and Evans 12 . Yates 13 , Iversen and 
Bernstein 14 , and Rosso w 15 . Although an inviscid model 
cannot describe the wake aging including its diffusion, it 
is still capable of producing the wake shape and its dy- 
namics. The mathematical modeling used in the above 
references were based on the use of the point vortex 
method to compute the motion of a finite number of point 
vortices which model the vorticity behind a wing. The 
first three-dimensional inviscid model was introduced by 
Kandil, et. al 16 , where the nonlinear vortex-lattice, which 
was also developed by Kandil, et. al 17 , was used to com- 
pute the interference flow between wings and the vortex 
wake hazardous effects. 

Viscous modeling of trailing vortices was first in- 
troduced by McCormick, et. al 1 *. Viscous interactions 
in vortex wakes and the effects of background turbu- 
lence, wind shear and ground on two-dimensional vor- 
tex pairs were presented by Bilanin, et. al 19, ®. Liu and 
his co-workers 21 ' 24 studied the interaction, merging and 
decay of vortices in two-dimensional space and of three- 
dimensional vortex filaments. To estimate the effects of 
density stratification, turbulence and Reynolds number on 
vortex wakes, an approximate model was recently devel- 
oped by Greene 25 . Lata on, Greene and his co-workers 26 
presented selected results of aircraft vortices which in- 
clude a juncture vortex, a lifting-wing vortex and a wake 
vortex. 

It is rrmrlmtftH from the above brief literature sur- 
vey that the problem contains several vortex flow regions 
along with several critical parameters influencing the vor- 
tex wake flows. The vortex wake flow of interacting-wing 
problems can be divided into three regions. The first re- 
gion includes the leading wing and its near-vortex-wake 
flow. The second region includes the interacting vortices 
in the wake. The third region includes the far-wake flow 
along with the trailing- wing flow. In Ref. 27, Kandil, 
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Wong and Liu presented results of the near-vortex-wake 
flow for a large-aspect-ratio rectangular wing using the 
thin-layer, Reynolds averaged, NS equations on a C-O 
grid. Next, a small aspect-ratio rectangular wing was in- 
troduced in the near-wake of the leading wing and the 
problem was recomputed to demonstrate the hazardous 
effects of the vortex wake flow. Two interference cases 
were considered. The first was called the “along-track 
penetration through vortex center” and the second was 
called the “along-track penetration between vortices”, see 
Fig. 1. It was demonstrated that in the first case, the 
trailing wing experienced large rolling moments and in 
the second case the trailing wing experienced loss of lift. 

In the present paper, the near-vortex-wake flow is 
recomputed using a finer grid than the grid used in Ref. 
27. Moreover, the thin-layer and full NS equations are 
used along with an upwind scheme, wherein the flux 
limiter of the solver is turned on and off. The goal of this 
part of the paper is to accurately compute the vortex wake 
flow and to study the effect of the numerical dissipation 
of the flux limiter on the computed results. In the second 
part of the paper, alleviation of the strength of the tip 
vortex without reducing the lift force is demonstrated by 
using directed flow injection from a slot along a certain 
length of the wing tip. 

FORMULATION 


Full and Thin-Layer Navier-Stokes Equations 

The conservative form of the dimensionless, unsteady, 
compressible, full Navier-Stokes equations in terms of 
time-independent, body-conformed coordinates p, p and 
P is given by 


dQ 

dE m 

0(E,-), ft , o 


dT + 

d£ m 

- —g^r 1 = 0; m=l-3,s = l- 3 

(1) 

where 






s m =£ m (xi,x 2 ,X3) 

(2) 


Q 

q l, lt 

= j = jb.PUj.puj.puj.pe] 

(3) 


E m = inviscid flux 

= jbUm.pUiUm + <9i£ m p,pU 2 U m 
+ ^£ m Pi£U3U m + d3£ m p,(/>e + p)U m ] 1 (4) 


(Ev) % = viscous and heat — conduction flux in 
direction 


dk^Krica - qk)]‘; k=l — 3, n=l — 3 


Urn = (6) 

The three momentum elements of Eq. (5) are given by 

W* £ Tte^*?** - !*■«*) fjs 

; j = 1-3 (7) 

The last element of Eq. (5) is given by 

V( u.% - <*) ■ ^ [( w 

1 d(a 2 )l 

(8) 

The single thin-layer approximations of the full Navier- 
Stokes equations demand that we only keep the deriva- 
tives in the normal direction to the body, p, in the viscous 
and heat flux terms in Eqs. (1), (7) and (8). Thus, we let 

s = 2 for the term in Eq. (1) and s = 2 and n = 2 
in Eqs. (7) and (8). These equations reduce to 


where 


dQ c?E» d(E,) a 

dt dp* dp ~ 



^ 2 (Uprk P -qk) = ^|v-W 

+*[^( U ? + U J + U 3 2 ) 

(T-l)P, dP J J 


(9) 

( 10 ) 


( 11 ) 


= d> eotf, i> = w = a p ^u p (12) 

In Eqs. (1M12), the dimensionless variables are refer- 
enced to their appropriate frees tream values. The dimen- 
sionless density p, cartesian velocity components ui, u 2 
and U), total energy per unit mass, e, dynamic viscosity, p 
and speed of sound, a, are defined as the ratio of the cor- 
responding physical quantities to those of the freestream; 
namely. P<*. Poe«2e. p<x and 0*; respectively. The 
pressure, p, is non -dime nsionalized by p^ a^, and is re- 
lated to the total energy far a perfect gas by the equation 
of state 


(5) 


P = (7 ~ - ^UjUj j ; j = 1 - 3 (13) 



where 7 is the ratio of specific heats and its value is 1.4. 
The viscosity, n , is calculated from the Sutherland law 

H = T 3/2 (y~“) t c = 0.4317 (14) 

where T is the temperature which is non-dimensionalized 
by T^. The Prandtl number, P r is fixed at 0.72. The 
Reynolds number is defined as R, = pocUocL/ l*<x nnd 
the characteristic length, L, is chosen as the root chord 
of the wing. 

In Eqs. (1)-(13), the indicial notation is used for 
convenience. The subscripts k, n, p and j are summation 
indices, the superscript or subscript s is a summation 
index and the superscript or subscript m is a free index. 
The partial derivative is referred to by d±. 

Turbulent Flow 

For the turbulent flow, the Navier-Stokes equations 
are transformed to the Reynolds-averaged equations by 
replacing the coefficient of molecular viscosity, /r, and 
the coefficient of thermal conductivity k with 


^ - ^(1 + l*t/p) 

(15) 


(16) 


where /r e is the effective viscosity, k« the effective ther- 
mal conductivity, ji, the turbulent viscosity, P, the lami- 
nar Prandtl number, P n the turbulent Prandtl number and 
C p the specific heat under constant pressure. The turbu- 
lent viscosity is obtained by using the two-layer al- 
gebraic eddy viscosity model which was first developed 
by Cebeci 2 ® for the boundary-layer equations and modi- 
fied later by Baldwin and Lomax 29 for the Navier-Stokes 
equations. The details of the turbulent model is given in 
Ref. (30) by Wong, Kandil and Liu. 

Boundary and Initial Conditions 

Boundary conditions are explicitly implemented. 
They include the inflow-outflow conditions and the 
solid-boundary conditions. At the plane of geometric 
symmetry, periodic conditions are used. The inflow- 
outflow boundary conditions are implemented using the 
one-dimensional Riemann-invariant conditions normal to 
these boundaries. On the solid boundaries, the no-slip 
and no-penetration conditions are enforced; Ui = U2 « uj 
= 0, and the normal pressure gradient is set equal to zero. 
For the temperature, the adiabatic boundary condition is 
enforced at the solid boundaries. 

The initial conditions correspond to the frees tream 
conditions with ut = U2 = uj = 0 on the solid boundaries. 


The freestream conditions are given by 

Poc = ^oc = Toe = 3> 
u joe = MocCoea, 

ujoc = MoeSina, (17) 

Ujoc = 0, 

- 1 / 1 , Ml 

Poc - 1/7>«0C - J) + 2 

where a is the angle of attack. 

COMPUTATIONAL SCHEME 

The implicit, upwind, flux-difference splitting, finite- 
volume scheme is used to solve the unsteady, com- 
pressible, thin-layer and Full Reynolds-averaged, Navier- 
Stokes equations. The scheme uses the flux-difference 
splitting scheme of Roe which is based on the solution of 
the approximate Riemann problem. In the Roe scheme, 
the inviscid flux difference at the interface of computa- 
tional cells is split into two parts; left and right flux dif- 
ferences. The splitting is accomplished according to die 
si gns of the eigenvalues of the Roe averaged-Jacobian 
matrix of the inviscid fluxes at the cell interface. The 
smooth flux limiter is used to eliminate oscillations at 
locations of large flow gradients. The viscous-and heat- 
flux terms are linearized and the cross -derivative terms 
are eliminated in the implicit operator. The viscous terms 
are differenced using a second-order accurate central dif- 
ferencing. The resulting difference equation is approx- 
imately factored and is solved in three sweeps in the 
^ and £ 3 directions. The computational scheme is 
coded in the computer program “FTNS-3D” 31 . 

COMPUTATIONAL RESULTS AND DISCUSSIONS 

We consider a rectangular wing of aspect ratio, AR 
= 5.9 and a NACA 0012 chordal section, at an angle of 
attack of o=4.64°. The flow Reynolds number, R« is 
3.2x10 s (based on the root-chord length) and the Mach 
number, M« is taken as 0.3 for low speed flows. The 
computational domain consists of a hemispherical bound- 
ary which is followed by a cylindrical boundary. The 
hemisphere center is located at the intersection of the 
wing root-chord and its trailing edge. Its radius is 15 
chords and the cylindrical boundary extends 14 chords in 
the downstream direction beyond the trailing edge. A C- 
O grid of 231x65x65 grid points in the stream wise wrap- 
around, spanwise and normal directions, respectively, is 
generated by using transfinite interpolation. The min- 
imum grid spacings normal to the wing surface at the 
leading-edge and tip regions are chosen to be 10" . The 
grid points close to the trailing edge are clustered up to 5 
chords in the downstream direction. The reason is that for 
the present work our interest is focused on the resolution 
of the near-vortex-wake flow. Figure 2 shows a typical 
grid for the wing and its near wake. The flow is consid- 
ered to be laminar until the chord station of 0.05 from 
the wing i^ing edge. Thereafter, the turbulent model is 



turned on to mimic the transition to turbulent flow. This 
location of the transitional flow has been experime ntall y 
determined by Yip and Shubert 32 . 

Thin-Layer and Full Navier-Stokes Solutions 

The flow around the rectangular wing is solved us- 
ing the thin-layer NS equations and the full NS equa- 
tions. The computations are carried out three times using 
the same grid described above and the implicit upwind 
scheme. In the first and second times, the thin-layer NS 
equations are solved using the implicit upwind scheme 
with and without the smooth flux limiter. In the third 
time, the full NS equations are solved using the implicit 
upwind scheme without the smooth flux limiter. Next, 
we compare the results of these three computions with 
each other and with the experimental data of Ref. 32. 

Figure 3 shows a comparison of the history of the 
logarithmic residual error and the lift coefficient of these 
computations. The thin-layer computation with a flux 
limiter (TL-limiter) shows a drop of the residual error of 
about five orders of magnitude after 8000 iteration steps. 
The thin-layer computation without a flux limiter (TL-no 
limiter) shows a drop of the residual error of the same 
order of magnitude as that of the TL-limiter computation 
after 10,800 iteration steps. The full Navier-Stokes com- 
putation without a flux limiter (FNS-no limiter) shows 
a drop of the same order of magnitude as that of the 
TL-limiter computation after 8,900 iteration steps. The 
lift coefficient of the three computations is 0.34090 (TL- 
limiter), 0.35354 (TL-no limiter) and 0.35250 (FNS-no 
limiter). In the next table, we show comparison of the 
three computations and the experimental data for the nor- 
mal force coefficient, Cn, the pitching-moment coefficient 
about the leading edge. Cm, and the drag coefficient, Co. 
The results of the case of next section for the tip-jet con- 
trol are also included in the table. The results show that 



C N 

Cm 

Cd 

TL-limiter 

0.34128 

-0.07581 

0.01857 

TL-no limiter 

0.35374 

-0.08161 

0.01682 

FNS-no limiter 

0.35272 

-0.08137 

0.01702 

Experiment 32 

0.35 

-0.0825 

N/A 

FNS-no limiter 
(tip jet control) 

0.3682 

-0.08667 

0.01732 


the Cn of the FNS-no limiter computation is the closest 
to the experimental C N with an error of +0.78%. The 
corresponding error for the TL-no limiter computation is 
+1.07% and for the TL-limiter computation is -2.49%. 
The error in the computed Cm in comparison with the 
experimental value is —1.37% for the NS-no limiter com- 
putation, -1.08% for the TL-no limiter computation and 
-8.11% for the TL-limiter computation. 

Figure 4 shows a comparison of the computed surface- 
pressure-coefficient distribution in the chordwise direction 
at different spanwise stations with the same experimental 


data. The present computed results are in good agreement 
with those of the experimental data with the exception of 
the peak suction pressure at the wing leading edge. This 
is attributed to the simple modeling of the transitional 
flow at this location, wherein a simple algebraic turbulent 
model is turned on. At the spanwise station of 0.9883 in 
the tip region, it is noticed that the computed Cp using 
the TL-no limiter and FNS-no limiter computations are 
in excellent agreement with the experimental data than 
that of the TL-limiter computation, particulary in the 
chordwise direction as of X/C=0.3. 

Figure 5 shows a comparison of the spanwise varia- 
tion of the total pressure-loss contours at chord-stations 
range of X/C=0.9-5.0 covering the evolution of vortex 
wake. The range of the total pressure-loss contours shows 
that the tip vortex reaches its maximum strength very 
close to the downstream side of the trailing edge. There- 
after, the tip- vortex core expands due to the viscous diffu- 
sion and moves in boards while growing in size. It is also 
observed that the wake thickness increases and moves up- 
wards. The TL-no limiter results of the vortex core and 
wake shear layer show less viscous diffusion and mere 
inboard motion than those of the TL-limiter results. The 
FNS-no limiter results show a little better resolution of 
the vortex core and wake shear layer than that of the TL- 
no limiter results. However, the inboard motion of the 
vortex core and wake shear layer is the same as that of 
the TL-no limiter results. Hence, it is concluded that the 
flux limiter is responsible for the vortex-core and wake 
shear-layer diffusions and their small inboard motions. In 
Fig. 6, the particle traces for the wing tip flow and the 
vortex wake flow are shown for the TL-limiter and FNS- 
no limiter computations. Again, the TL-limiter results 
show the larger diffusion of the tip vortex than that of the 
FNS-no limiter results. 

Active Injection Control of Tip Vortex 

The hazardous effects of the tip vortex on trailing 
aircraft is very significant as it has been demonstrated in 
Ref. 27. Several ideas exist for passive and active control 
of the tip vortex in order to alleviate its hararrimic effects. 
The challenging issue here is to diffuse the tip vortex 
without decreasing the lifting force of the leading aircraft. 
Our first a tt e mp t to accomplish these objectives is to use 
active injection from a slot along the wing tip. The slot 
is located on the wing-tip surface and extends from the 
chord station X/C-0.2 to the chord station X/C=0.6. The 
direction of the jet blowing is downstream at 45® which 
is measured from the x-axis. A constant mass-flow-rate 
coefficient of 0.015 is used for the jet blowing through 
the wing-tip slot. The solution of this case is obtained 
by using the FNS-no limiter computation. Two types 
of initial conditions have been used. In the first, the 
control starts from the uniform flow conditions while in 
the second, the control starts from the FNS-no limiter 
solution of the preceding section. The final steady flow 


solution for both cases showed the same results. Here, 
only the results of the first case are shown. 

Figure 7 shows the history of the residual error and lift 
coefficient for 16,000 iteration steps. The spikes shown 
in the residual error graph are due to the step increases 
of the mass-flow rate from 0.005 to 0.01 and finally to 

0.015. The residual error drops three orders of magnitude 
and the lift coefficient shows a value of 0.3680. The lift 
coefficient is 4.4% higher than that of the FNS-no limiter 
computation without injection. The values of Cn, Cm 
and Cd are included in the table of the preceding section. 
The Cd coefficient shows 1.76% increase than that of the 
FNS-no limiter computation without injection. In Fig. 8, 
the pressure coefficient is shown at the spanwise stations 
near the wing tip for comparison with the experimental 
data of the case without injection. It is observed that 
only the tip-section C p near the trailing edge shows a 
slight decrease from that of the experimental data. 

Figure 9 shows the spanwise variation of the total 
pressure-loss contours at chord stations from X/C=0.95 
to X/C=5.0. A comparison of these results with those 
of Fig. 5 reveals the excellent effect of the tip-injection 
control in dissipating the tip vortex. This is also seen in 
Fig. 10 of the particle traces in comparison with those 
of Fig. 6. 

Figure 11 shows a blow-up of the total pressure- 
loss contours at the wing tip covering chord stations 
from X/C=0.2 to X/C=0.9. It is observed that the tip 
vortex existing at X/C=0.2 (at the starting location of the 
slot) is quickly diffused by the blown jet. However, by 
the end of the slot length at X/C=0.6, the tip vortex is 
recovering but with small strength. It is concluded that 
tip-injection control is effective in diffusing the tip vortex 
and moreover it is very important to study the effect of the 
injection-slot length in the future. Other parameters need 
to be included in the future study such as the direction of 
injection, the mass-flow rate and the width of injection. 

CONCLUDING REMARKS 

The thin-layer and full NS equtions have been used 
to accurately solve for the near-vortex-wake flow around 
a large-aspect-ratio rectangular wing. The flow has been 
computed by using the thin-layer equations with and with- 
out flux limiters and the NS equations without flux lim- 
iters, all on the same grid. The results show the substan- 
tial effects of the flux limiter in diffusing the tip vortex 
and the trailing-edge shear layer and in their small in- 
board motions. The full Navier-Stokes computation with- 
out flux limiters show the best results in comparison with 
the experimental data. Also, it has been demonstrated that 
tip-injection control is substantially effective in diffusing 
the tip vortex without reduction of the wing lift coef- 
ficient Although the drag coefficient slightly increased 
the lift to drag ratio increased. Currently, work is under- 
way to conduct an extensive parametric study on the tip- 
injection control and to study its effect on a trailing wing. 
Passive-control methods are also being investigated. 
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Figure 1. Types of interference of Trailing Voitex Wake; (1) cross-track penetration, (2) 
along-track penetration between dp vortices, and (3) along-track penetration 
through dp vortex center. 



Figure 2. Typical computational grid for the wing 
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Figure 3. Iteration history of residual error and lift coefficient 
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Figure 4. Chordwise surface-pressure distribution at spanwise stations for an isolated 
wing and comparison with experimental data. 
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Figure 7. Iteration history of residual error and lift coefficient for the wing with tip- 
injection control (F.N.S. no limi ter) 



Figure 8. Chord wise surface-pressure distribution at span wise stations for an isolated 
wing with tip-injection control (F.N.S. no limiter). 



Figure 10. Particle traces for the wing tip with tip-injection control and its vortex wake 
(F.N.S.-no limiter). 
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Figure 11. Span wise variation of total pressure-loss contours at chord wise stations on the 
wing with tip-injection control (F.N.S.-no limiter). 
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ABSTRACT 

The effectiveness of certain active-control methods 
for asymmetric flows around circular cones is investigated 
by using computational solution of the unsteady, com- 
pressible full Navier-Stokes equations. Two main meth- 
ods of active control which include flow injection and 
surface heating are used. For the flow-injection-control 
method, flow injection is used either in the normal di- 
rection to the surface or in the tangential direction to 
the surface. For the surface-heating-control method, the 
temperature of the cone surface is increased. The effec- 
tiveness of a hybrid method of flow control which com- 
bines normal injection with surface heating has also been 
studied. The Navier-Stokes equations, subjected to vari- 
ous surface boundary conditions, are solved by using an 
implicit, upwind, flux-difference splitting, finite-volume 
scheme for locally-conical flow solutions. 


INTRODUCTION 


The problems of prediction, analysis and control of 
asymmetric vortical flows around slender pointed bodies 
are of vital importance to the dynamic stability and con- 
trollability of missiles and fighter aircraft. The onset of 
flow asymmetry occurs when the relative incidence (ratio 
of angle of attack to nose semi-apex angle) of pointed 
forebodies exceeds certain critical values. At these crit- 
ical values of relative incidence, flow asymmetry devel- 
ops due to natural and/or forced disturbances. The origin 
of natural disturbances may be a transient side slip, an 
acoustic disturbance, or similar disturbance of short du- 
ration. The origin of forced disturbances is geometric 
perturbations due to imperfections in the nose geometric 
symmetry or similar disturbances of permanent nature. 
In addition to the relative incidence as one of the influ- 
ential parameters for the onset of flow asymmetry, the 
freestream Mach number, Reynolds number and shape of 
the body-cross sectional area are also important param- 
eters. 

In several recent papers by Kandil et aL M , the un- 
steady, thin-layer, compressible Navier-Stokes equations 
have been used to simulate steady and unsteady, asym- 
metric vortex flows, including their passive control. 
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around cones with different cross-sectional shapes. The 
emphasis of these papers was extensive computational 
studies of the parameters which influence the asymmet- 
ric flow phenomenon and its passive control. Since the 
computational cost associated with the solution of three- 
dimensional-flow problems at reasonable flow resolution 
is very expensive, all the computational solutions were 
obtained using a locally-conical flow assumption. Such 
an assumption reduces the problem solution to that on 
two conical planes, which are in close proximity of each 
other, and hence it reduces the computational cost by an 
order of magnitude. Moreover, such solutions still pro- 
vide extensive understanding of the flow physics since 
one can use very fine grids for reasonable flow resolu- 
tion. These studies showed that asymmetric flow solu- 
tions were unique irrespective of the type of flow distur- 
bance; a random disturbance in the form of a machine 
round-off error or a controlled disturbance in the form 
of a short-duration side-slip disturbance. Unsteady asym- 
metric flow solutions with perfectly periodic vortex shed- 
ding were successfully simulated, and the solutions were 
unique irrespective of the computational scheme used. It 
has also been shown that as the Mach number was in- 
creased, the flow asymmetry was decreased, and as the 
Reynolds number was increased the flow asymmetry was 
increased. Moreover, the cross-sectional shape of the 
cone has been shown to be a very influential parameter 
on the flow asymmetry. Circular sections produced very 
strong flow asymmetry and diamond sections produced 
rclatively-weaker flow asymmetry. 

In a later paper, by Kandil et aL s , the full Navier- 
Stokes solutions were compared with the thin-layer 
Navier-Stokes solutions. It was shown that the full 
Navier-Stokes solutions produced thicker free-shear lay- 
ers and more vortex-core resolution as compared with 
those of the thin-layer Navier-Stokes equations. In Ref. 

5, a few tentative three-dimensional flow solutions were 
also presented. 

Substantial research efforts have recently been de- 
voted for eliminating or alleviating flow asymmetry 
and its corresponding side force. In the experimen- 
tal area, several passive-control methods 64 and active- 
control methods*' 13 have been investigated. Computa- 
tional simulations have also been used to investigate the 
effectiveness of several passive-control methods 13 and 
active-control methods 12 * w * 15 . Various methods of pas- 
sive control were demonstrated in the above references 
which include the use of vertical fins along the leeward 
plane of geometric symmetry, thin and thick side strokes 


with different orientations, and rotatable forebody tips 
which have variable cross section (from a circular shape 
at its base to an elliptic shape at its tip). It was shown by 
Kandil et al. 4 that side-strakes control is more practical 
than the vertical-fin control since the former was more ef- 
fective over a wide range of angle of attack than the latter. 
Moreover, side-strake control provided an additional lift- 
ing force. However, the effectiveness of the side-strake 
control terminates at very high angles of attack for the 
considered strake geometry and flow conditions. 

Various active-control methods have been used which 
include forebody blowing and movable forebody strakes. 
The forebody blowing methods include forward blowing, 
normal blowing, aft blowing and tangential blowing. The 
main concept of forebody blowing is to control flow 
separation on the forebody and to c re a te yawing forces 
and moments which can be utilized in controlling the 
body. 

In the present paper, we investigate the effectiveness 
of two main methods of active control which include flow 
injection and surface heating. The study of flow-injection 
control covers normal and tangetial injection. Moreover, 
a hybrid method of flow control which combines surface 
heating and normal injection methods is also investigated. 
At this stage of research, the flow control is aiming at 
eiftM-ir rendering the vortical asymmetric flow symmetric 
or rendering the surface-pressure distribution symmetric. 
Computational solution of the unsteady, compressible, 
full Navier-Stokes equations is used for the present work 
with the exception of the tangential injection. For the 
tangential injection, the thin-layer NS equations are used. 
The computational applications are focused on circular- 
section cones, and locally -conical flow assumptions are 
used to substantially reduce the computational cost 

FORMULATION AND COMPUTATIONAL SCHEME 


Full Navier-Stokes Equations 

The conservative form of the dimensionless, unsteady, 
compressible, full Navier-Stokes e qu atio ns in terms of 
time-independent, body-conformed coordinates (' , and 
£ 3 is given by 
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The three momentum elements of Eq. (5) are given by 

The last element of Eq. (5) is given by 

;p ■ 1 - 3(81 

The reference parameters for the dimensionless form 
of the equations are L, a^.L/a, x ,p x and fa the 
length, velocity, time, density and molecular viscosity, 
respectively. The Reynolds number is defined as Re * 
PoeVoeL/^oe, where the characteristic length, L, is the 
body length. The pressure, p, is related to the total en- 
ergy per unit mass and density by the gas equation 

P = (7~ l)p|«~ ^( u i + u i + u j)| (9) 
The viscosity is calculated from the Sutherland law 

A* = tS/1 . c = 04317 < 10 > 

and the Prandtl number P r m 0.72. In Eqs. (1)-(10). the 
indidal notation is used for convenience. 

Boundary and Initial Conditions 

The boundary and «"«*«> conditions vary according 
to the problem under consideration. The boundary con- 
ditions are explicitly satisfied. In general, they include 
inflow-outflow conditions and solid-boundary conditions. 
For problems of flow asymmetry, where the flow is solved 
throughout the whole computational domain, periodic 
boundary rnmtitinns ^ used at the plane of geometric 
symmetry of the problem. 

For the asymmetric flow problems around slender 
bod i es and for supersonic inflow-outflow boundary, the 
Riemann-invariant boundary conditions are used. They 
require that the inflow variables be at the frees tream con- 
ditions, and the conical shock enclosing the body be 


captured as part of the solution. For supersonic out- 
flow boundary, the Riemann-invariant boundary condi- 
tions require that all flow variables be extrapolated from 
the interior cells. On the solid boundary, without injec- 
tion or heating, the no-slip and no- penetration conditions 
are enforced. Moreover, the zero normal-pressure gradi- 
ent and adiabatic boundary conditions are enforced. For 
the present active control problems, the mass-flow rate 
is specified at the body surface for the normal injection 
control and the temperature distribution is specified at 
the surface for the heating control. For the tangential 
flow injection the mass flow rate and velocity profile are 
specified at the lip exit (shown in Fig. 10 of next section). 

The initial conditions correspond to the uniform flow 
conditions with ui = U2 = U3 = 0 on the solid boundary. 
These conditions are used to obtain the asymmetric flow 
solution. Next, the flow control conditions are enforced 
and the previously obtained asymmetric solution is used 
for the initial conditions of the active control problem. 

Computational Scheme 

The implicit, upwind, flux-difference splitting, finite- 
volume scheme is used to solve the unsteady, compress- 
ible, full Navier-Stokes equations. The scheme uses the 
flux-difference splitting scheme of Roe which is based on 
the solution of the approximate Riemann problem. In the 
Roe scheme, the in viscid flux difference at the interface of 
computational cells is split into two parts; left and right 
flux differences. The splitting is accomplished accord- 
ing to the signs of the eigenvalues of the Roe averaged- 
Jacobian matrix of the in viscid fluxes at the cell interface. 
The smooth flux limiter is used to eliminate oscillations at 
locations of large flow gradients. The viscous-and heat- 
flux terms are linearized and the cross-derivative terms 
are eliminated in the implicit operator. The viscous terms 
are differenced using a second-order accurate central dif- 
ferencing. The resulting difference equation is approxi- 
mately factored and is solved in three sweeps in the 
£ 2 , and £ 3 directions. The computational scheme is coded 
in the computer program ‘TTNS3D". 

For the locally-conical flow solutions, an axial station 
of xi =« 1.0 is selected and the components of the flowfield 
vector are forced to be equal between this axial station and 
another axial station in close proximity to xi = 1.0. This 
ensures that the flow variables are locally independent of 
the axial direction at xi » 1.0 (Kandil et al 5 ). 

COMPUTATIONAL APPLICATIONS AND DISCUSSION 

For all the computational applications shown in this 
section, a 5°-semiapex angle circular cone at 1.8 Mach 
number and 10 5 Reynolds number is considered. These 
flow conditions were considered earlier (Kandil et al. 4, 5 ) 
for the same cone for asymmetric-flow prediction and 
passive-control methods. The grid is 241x81x2 points in 
the wrap around, normal directions and axial direction, 
respectively. The grid is generated using a modified 


Joukowski transformation with a minimum grid size of 
ICT 4 in the £ 2 direction at the body surface. For the 
tangential flow injection a multi-block grid has been used 
and it is explained later on. 

Asymmetric Flow, a = 20° 

Figure 1 shows the locally-conical flow solution with- 
out any control around the cone at an angle of attack 
of 20°. The figure shows the total-pressure-loss (TPL) 
contours and the surface-pressure (SP) coefficient. The 
surface pressure is presented versus the angle 6, which 
is measured in the clockwise direction from the leeward 
plane of geometric symmetry. This stable asymmetric 
flow solution is obtained after 6,000 iteration steps. The 
source of flow disturbance at these critical flow condi- 
tions is the truncation error. As it has been shown earlier 
(Kandil et al. 1 ), this solution is unique irrespective of the 
type of source of the flow disturbance. 

Normal Injection Control, 9 = - 22J° — ► - 67-5°, 
m = constant = 0.03, a * 20° 

Next, a constant mass-flow injection of rti * 0.03 is 
applied normal to the cone surface. The circumferential 
range of injection ports extends over $ = -22.5® — ► -67.5°. 
The solution of the previous asymmetric flow case is 
used as initial conditions along with modified surface 
boundary conditions. Figure 2 shows the TPL contours 
and the SP coefficient from the solution of the controlled 
flow. Although the vortical flow is still asymmetric, as 
it is seen from the TPL contours, the SP curve shows 
a symmetric distribution resulting in a zero side force. 
The TPL contours show that a primary vortex of smaller 
height (in comparison to the case of Fig. 1) is still existing 
on the right side and it is connected through a free-shear 
layer with the left side of the body surface. The injection 
flow on the left side of the body decreases the suction 
pressure on that side making its distribution equal to that 
on the right side. 

Normal Injection Control, 0 = 0° — - 67.5°, 
ib s constant = 0.03, a = 20° 

In this case, the constant mass-flow injection of rti 
= 0.03 is kept fixed while the circumferential range of 
injection ports is extended to cover the range of 6 * 0° 
— ► -67.5°. Figure 3 shows the TPL contours and SP 
coefficient of the solution of the controlled flow. The 
vortical flow is still asymmetric but the SP curve stows a 
symmetric distribution, which results in a zero side force. 
It is also noticed that the primary vortex on the right side 
moves to a larger height (in comparison with the cases of 
Figs. 1 and 2) and the left-side free-shear layer also moves 
to the same height level. The SP coefficient curve stows 
less negative pressure coefficient on the leeward side 
in comparison with the case of Fig. 2, which produces 
smaller lifting force. 


Normal Injection Control, 9 = - <7i* — + 47.5®, 
Variable Mam-Flow Rate, in, * 0.03, q ■ 20° 

In this case, the circumferential range of injection 
ports is extended to cover a symmetric range of 9 ■ - 
67.5® — +67.5®. Moreover, the injected mass-flow rate 
is made proportional to the difference in the surface pres- 
sure between corresponding points on the right and left 
side of the body (within 9 = - 67.5°— ► +67.5°). The maxi- 
mum injected mass-flow rate, rtw. which corresponds to 
the maximum pressure difference is restricted to ihn*, = 
0.03. Practically, this control method can be achieved by 
sensing the pressure difference between the right and left 
ports and using it as a feedback control in order to inject 
a mass-flow rate which is proportional to that pressure 
difference. Figure 4 shows the results of this controlled 
flow. The history of the residual error versus the number 
of iterations shows the asymmetric solution up to 6,000 
iterations and the evolution of the symmetric controlled- 
flow solution up to 7,200 iterations. The corresponding 
curve of the history of the side force shows that the side 
force becomes zero at the end of 7,200 iterations. Snap- 
shots of the evolution of the symmetric controlled flow 
solution, in terms of the TPL contours, are given at the 
6,200; 6,400; 6,600; 6,800; 7,000 and 7,200 iteration step. 
The vortical flow is rendered symmetric using the present 
method and the SP coefficient curve shows a symmetric 
distribution. 

Surface Beating Control, 9 ■ • — -180®, 

Ti * 7 Toe, a * 20° 

For the heating control, the left surface temperature 
of the cone, in the circumferential range of $ = 0® — ► 
-180®, has been changed from the adiabatic condition to 
a constant temperature condition of surface temperature, 
T, = 7 Toe, where Toe is the freestream temperature. 
Starting from the stable asymmetric flow solution at the 
12,000 iteration step, the left surface temperature is raised 
to T, = 5Toe and then it is raised further to T, ■ 7 Toe. 
Figure 5 shows the results of this case which include 
the history of the residual error and the corresponding 
history of the side force along with the TPL contours 
after 24,000 iterations. The side-force curve shows that 
the side force decreases as the temperature increases until 
it vanishes for T,« 7 T*. Although the side force is zero, 
the vortical flow is still asymmetric. This case shows the 
effect of increasing the gas viscosity with heating in order 
to equalize the surface pressure distribution on the right 
and left sides of the body. This is the first time, we know 
of. that such an active control method has been applied 
to control asymmetric flows. 


Iqjectioa Flow Control, Variable Maaa-Flow 
Rate, a « 20° - 30® 

In this case, the effectiveness of the normal injection 
control of the case of Fig. 4 (rti_. a 0.03, a ■ 20®, 9 = 

- 67.5* — ► 67.5*) is studied for higher angles of attack. 
Figure 6 shows the results of this study which includes 
the history of the residual error versus the number of 
iterations, the corresponding history of the side force and 
the TPL contours along with the SP-coefficient curves for 
a = 22®, 24®, 26®, 28® and 30®. The residual error and the 
side- force curves show the responses for the following 
history: Fust, no control is applied until the 12,000 
iteration step; second, injection is applied at o = 20® 
for 2,000 iterations; third, the angle of attack is increased 
2* each 500 iterations until a =* 28°; fourth, at o « 30° 
injection is continued for 2,000 iterations. The side-force 
curve shows zero value in the angle of attack range of a = 
20® - 28* (corresponding to the iterations range of 14,000 

- 16,000). When the angle of attack reaches 30®, the 
normal injection control becomes incapable of achieving 
flow symmetry. It is interesting to study the snapshots of 
the TPL contours at these angles of attack. It is noticed 
that the primary vortices increase in size in the normal 
direction and their inner boundaries approach each other. 
At a * 28®, the inner boundaries of the primary vortices 
touch each ocher, and thereafter at a « 30® the primary 
vortices become asymmetric. The asymmetric response 
of the primary vortices at a * 30® is believed to be due 
to a strong instability arising due to the strong interaction 
of the inner shear layers of the primary vortices. 

At a * 30®, the injected mass-flow rate in the normal 

direction to the body surface is increased to rh * 0.05, 

0.06 and 0.07 in order to recover the flow asymmetry. 
Figure 7 shows the results of these tests. Both the TPL 
contours and the SP-coefficient curves show that the flow 
is asymmetric and the side force does not vanish. It is 
seen that the flow asymmetry changes from the right to the 
left sides and so does the side force. Therefore, it is ap- 
parent that normal-injection control loses its effectiveness 
at high angles of attack for the considered conditions. 

Hybrid H«ating-Iid«ctioa Control, a * 20° - 42® 

In this case, the surface-heating control is followed 
by normal-injection control. Figure 8 shows the results 
of this study. The control process is applied as follows; 
Fust, the whole surface is heated to T, = 5T« and the 
heating control started at a ■ 20®, where the number 
of iterations is 12,000. Next, the angle of attack is 
increased 4* until symmetric flow is recovered at the 
18,000 iterations. The angle of attack is then increased 
2* each 1000 iteration until the angle of attack reaches 
38*. The flow symmetry has been continuously recovered 
until the angle of attack of 36* where a slight flow 
asymmetry appears in the solution. At a - 38®, the flow 
asymmetry becomes stronger and the heating control loses 
its effectiveness for T, = 5 Toe- 


At o = 38°, normal-injection control is applied with 
iti — = 0.05 and 9 = - 67.5° — 67.5°. The angle of 
attack is then increased 2° each 3000 iterations. The 
history of the residual error and the corresponding zero 
side force is shown in Fig. 8. Figure 9 shows snapshots 
of the TPL contours and SP-coefficient curves at o = 38°, 
40° and 42°. The TPL contours show asymmetric flows 
but the SP-coefficient curves show symmetry resulting 
in zero side force. Thus, it is concluded that a hybrid 
control of heating and normal injection is very effective 
in eliminating the side force at very high angles of attack. 

Tangential-Injection Control, iti-., = 02 

In this case, the circular section of the lower portion 
of the cone has a raidus larger than that of the upper 
portion. The lower local radius is 1 .05 that of the upper 
one. Flow is injected from the lip of the lower portion in a 
direction tangent to the surface of the upper portion. The 
maximum mass flow rate, rtw, is 0.2 and the velocity 
profile at the lip exit is a parabolic one. Figure 10 shows 
the grid used for the thin-layer computations. Three 
grid blocks are used for the present computation. The 
first block is 121x19x2, the second is 121x63x2 and the 
third is 121x63x2 all in the circumferential, normal and 
axial directions, respectively. The tangential flow control 
starts from the 20° angle of attack. The computations are 
carried out until a symmetric flow solution is obtained. 
Thereafter, the angle of attack is increased to 30° with 
a step of 2°. Figure 10 also shows the total pressure- 
loss contours for the controled, symmetric flow solutions 
at a = 22°, 26° and 30°. It is noticed that as the 
angle of attack increases the inner shear layers of the 
controled symmetric primary vortices approach each other 
and hence flow asymmetry might develop again. 

CONCLUDING REMARKS 

The computational solution of the unsteady, com- 
pressible, full Navier-Stokes equations has been used to 
study the effectiveness of two active control methods for 
asymmetric flows around circular pointed cones. Locally- 
conical flow assumptions have been used in order to re- 
duce the computational cost of this study by an order of 
magnitude in comparison with that of three-dimensional 
flow solutions. The first active control method is applied 
by injecting air normal to the body surface from ports that 
are circumferentially distributed within a certain angle. 
The injected mass-flow rate is either constant or variable. 
The variable injected mass-flow rate is proportional to the 
surface-pressure difference between corresponding points 
on the right and left sides of the body surface. The results 
of this part of the study show that use of normal injection 
alone loses its effectiveness at very high angles of attack. 
The second active control method is applied by heating 
the body surface either partially or totally. The results of 
this part of the study show that surface-heating control is 
much more effective in comparison with normal-injection 


control at very high angles of attack. Next, a hybrid con- 
trol method which includes surface-heating control fol- 
lowed by normal-injection control is applied. The results 
of this part of the study shows the superior control ef- 
fectiveness of the hybrid method. Active control using 
tangential flow injection is very promising. Work is still 
underway to exploit its effectiveness. Simultaneous ap- 
plication of passive control and the present active control 
methods is also investigated. 
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coefficient, a = 30°, Moo = 1.8, R« = 10 s , variable 
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Figure 8. Hybrid heating- injection control (heating part); 

residual error and side force and TPL contours and 
SP coefficient; a = 20° - 38°, Moo = 1.8, R« = 10 5 , 
T. = 5 Too on whole surface. 
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Figure 9. Hybrid heating-injection control (injection part); 

TPL contours and SP coefficient; a = 38* - 42*! 
Mo,, = 1.8. R, = 10 5 , variable rh. iti^, = 0.05. 9 = 
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Figure 10. Typical multi -block grid and total pressure-loss contours for the 
tangential-flow injection control; a = 22°, 26°. 30°; Moo = 1.8, Re 
= 10 5 , parabolic velocity profile at lip exit, rh~.. = 0.2 
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ABSTRACT 

The unsteady, compressible, full Navier-Stokes (NS) 
equations are used to study the critical effects of the 
downstream boundary conditions on the supersonic vortex 
breakdown. The present study is applied to two super- 
sonic vortex breakdown cases. In the first case, quasi- 
axisymmetric supersonic swirling flow is considered in a 
configured circular duct, and in the second case, quasi- 
axisymmetric supersonic swirling jet, that is issued from 
a nozzle into a supersonic jet of lower Mach number, is 
considered. For the configured duct flow, four different 
types of downstream boundary conditions are used, and 
for the swirling jet flow from the nozzle, two types of 
downstream boundary conditions are used. The solutions 
are lime accurate which are obtained using an implicit, 
upwind, flux -difference splitting, finite-volume scheme. 

INTRODUCTION 

Numerical simulation of vortex breakdown has been 
focused in most of the existing literature on incompress- 
ible flows. Quasi-axisymmctric, vortex-breakdown sim- 
ulations of incompressible swirling flows have been pre- 
sented by Grabowski and Berger 1 , Hafez, et al. 2 , Salas 
and Kuru villa 3 , Menne 4 , Wu and Hwang 3 and Menne and 
Liu 6 . Three-dimensional, vortex-breakdown simulations 
of incompressible swirling flows have been presented by 
Spall, et al. 7 , Hanel*, and Krause 910 . Discussion of the ef- 
fects of side boundary conditions and downstream bound- 
ary conditions on vortex breakdown was presented by 
Krause in Ref. 10. 

Longitudinal vortex/lransverse shock-wave interac- 
tions are typical applications which appear in transonic 
and supersonic flows over a strake-wing configuration 
at moderate-high angles of attack, at a supersonic in- 
let injesting a vortex and inside a supersonic combustor 
where fuel is injected in a swirling jet to enhance fuel- 
air mixing 1113 . For the strake-wing configuration, vortex 
breakdown is undesirable since it results in the stall phe- 
nomenon, and hence its occurrence need to be delayed. 
On the other hand, vortex breakdown for the other two 
applications is desirable since it enhances mixing and sta- 
bility of the flame 14,13 , and hence its occurrence need to be 
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controlled for the optimum performance. Unfortunately, 
the literature lacks this type of analysis with the excep- 
tion of the preliminary work of Liu, Krause and Menne 16 , 
Copening and Anderson 17 , Delery, et al, 11 Kandil and 
Kandil 1 * and Meadows, Kumer and Hussaini 19 . 

The first time-accurate NS solution for a supersonic 
vortex breakdown was developed by the present au- 
thors in Ref. 20. We considered a supersonic quasi- 
axisymmetric vortex flow in a configured circular duct 
The time-accurate solution of the unsteady, compress- 
ible NS equations was obtained using an implicit, up- 
wind, flux-difference splitting finite-volume scheme. A 
shock wave has been generated near the duct inlet and 
unsteady vortex breakdown has been predicted behind 
the shock. The predicted flow was characterized by the 
evolution, convection and shedding of vortex breakdown 
bubbles. The Euler equations were also used to solve the 
same problem. The Euler solution showed larger size and 
number of vortex-breakdown bubbles in comparison with 
those of the NS solutions. The time-accurate solution was 
carried out for 3,200 times steps which are equivalent to 
a dimensionless time of 16. Only one value of Reynolds 
number of 10,000 was considered in Ref. 20. 

In a later paper 21 , we expanded our study of this flow 
using time-accurate computations of the NS equations 
with a fine grid in the shock-vortex interaction region 
and for long computational times. Several issues were 
addressed in that study. First, we showed the effect of 
Reynolds number on the temporal evolution and persis- 
tence of vortex-breakdown bubbles behind the shock. In 
that stage of computations, the conditions at the down- 
stream exit were obtained by extrapolating the compo- 
nents of the flowfield vector from the interior cell centers. 
Although the flow was supersonic over a large portion of 
the duct exit, subsonic flow existed over a small portion 
of the exit around the duct centerline. Therefore, selected 
flow cases were computed using a Riemann-invariant- 
type boundary conditions at subsonic points of the duct 
exit Finally, the effect of swirl ratio at the duct inlet 
was investigated. 

In the present paper, we address the problem of spec- 
ifying the downstream boundary conditions and their crit- 
ical effects on the supersonic vortex breakdown problem 
for internal and external flows. For this purpose, the un- 
steady, compressible, full NS equations are used along 
with an implicit, upwind, flux-difference splitting, finite- 
volume scheme for the time-accurate solutions. For the 


inicmal flow case, supersonic swirling flow in a config- 
ured duct is considered along with four types of down- 
stream boundary conditions. Keeping the duct geometry 
and the upstream flow conditions fixed, the exit boundary 
conditions are varied. The four exit boundary conditions 
include extrapolation of all the five variables from the 
interior cell centers, specifying the downstream pressure 
and extrapolating the other flow conditions from the in- 
terior cell centers, specifying the downstream pressure 
gradient and extrapolating the other four conditions from 
the interior cell centers, and using a disk of specified 
radius at the exit section. For the external flow case, 
supersonic swirling jet is issued from a nozzle into a su- 
personic non-swirling jet of Mach number lower than that 
of the swirling jet. Two types of downstream boundary 
conditions are considered. In the first type, extrapolation 
of all the five variables from the interior cell centers is 
used, while in the second type, the standard Riemann- 
invariant type boundary condition is used. All the results 
presented in this paper have been produced using our full 
NS solver which is known as the “FTNS3D” code. The 
CRAY-YMP computer of NASA Langley Research Cen- 
ter is used for the computations. 

HIGHLIGHTS OF THE FORMULATION 
AND COMPUTATIONAL SCHEME 

The conservative, unsteady, compressible, full Navier- 
Slokcs equations in terms of lime-independent, body- 
conformed coordinates £ l , ( 2 and £ 3 are used to solve the 
problem. The equations are given in Ref. 21 and hence 
they are not shown here. Along with these equations, 
boundary conditions are specified at the computational- 
domain inlet, side wall and downstream exit. The down- 
stream exit boundary conditions will be presented in the 
next section of the computational results. The initial con- 
ditions will also be presented in the next section. 

The computational scheme used to solve the unsteady, 
compressible full NS equations is an implicit, upwind, 
flux -difference splitting, finite-volume scheme. It em- 
ploys the flux -difference splitting scheme of Roe which is 
based on the solution of the approximate one-dimensional 
Riemann problem in each of the three directions. In the 
Roe scheme, the inviscid flux difference at the interface 
of a computational cell is split into left and right flux dif- 
ferences. The splitting is accomplished according to the 
signs of the eigenvalues of the Roe averaged-Jacobian 
matrix of the inviscid flux at the cell interface. The 
smooth limiter is used to eliminate oscillations in the 
shock region. The viscous and heat-flux terms are lin- 
earized and the cross-derivative terms of the viscous Ja- 
cobians are dropped in the implicit operator. These terms 
arc differenced using second-order spatially accurate cen- 
tral differencing. The resulting difference equation is ap- 
proximately factored and is solved in three sweeps in the 

and directions. The scheme is used for third- 
order spatial accuracy and first-order temporal accuracy. 


The scheme is coded in the computer program which is 
called “FTNS3D". 

The quasi-axisymmetric solutions are obtained using 
the three-dimensional code by forcing the flowfield vector 
to be equal on two axial planes, which are in close 
proximity of each other. 


COMPUTATIONAL RESULTS AND DISCUSSION 


I. Vortex Breakdown in Configured Circular Duct 


Figure 1 shows a configured circular duct which con- 
sists of a short, straight cylindrical part at the inlet which 
is followed by a short divergent cylindrical part until the 
axial length of 0.74. The divergence angle is 6°. The duct 
radius is then kept constant and a convergent-divergent 
nozzle with a throat radius of 0.95 is attached. The duct 
exit radius is 0.98 and its total length is 2.9. The diver- 
gent part of the duct ensures the stability of the formed 
shock in the inlet region. The configuration of the duct is 
intended to ensure that the supersonic inflow will become 
supersonic at the exit. As the computations will show, a 
small portion of the duct exit near its centerline becomes 
subsonic at certain times for the specified inflow condi- 
tions. This configured duct has also been used by Delery, 
ct al. 11 for their Euler equations computations of super- 
sonic vortex breakdown in an attempt to computationally 
model an experimental set up. 


The NS solver is used to compute all the following 
flow cases by using a grid of 221 x51 on two axial planes, 
where 221 points are in the axial direction and 51 points 
are in the radial direction. In the inlet region up to 
the 0.74 axial station, 100 grid points are used and the 
remaining 121 points are used in the remaining part of the 
duct. The grid is also clustered at the centerline (CL) and 
the wall. The minimum radial grid size at the CL is 0.002. 
The two axial planes are spaced circumferentially at a 
certain angle so that the aspect ratio of the minimum grid 
size will be less than 2. The upstream Mach number is 
kept at 1.75 and the Reynolds number is kept at 100,000. 
The initial profile for the tangential velocity is given by 



( 1 ) 


where Ucc = 1.74, r* = 0.2 and it, = 0.1. The max- 
imum jS-, swirl ratio /?, is at r * 0.224 and its value is 
kept at 6.32. The radial velocity, v, at the initial station 
is set equal to zero and the radial momentum equation is 
integrated to obtain the initial pressure profile. Final ly 
the density p is obtained from the definition of the speed 
of sound for the inlet flow. With these compatible set 
of profiles, the computations are carried out accurately in 
time with At = 0.0025. The wall boundary conditions fol- 
low the typical Navier-Stokes solid-boundary conditions. 
These computations have been carried out on the CRAY 
YMP of the NASA Langley Research Center. The CPU 
time is 40 ps/grid point/iteration for the NS calculation. 



Next, we present the results of the computational 
study of this case which covers four types of exit bound- 
ary conditions. 

1.1. Extrapolation from Interior Cell Centers 

Figure 2 shows snapshots of the streamlines and Mach 
contours for this case. The downstream boundary condi- 
tions are obtained by extrapolating all the flow variables 
from the cell centers at the exit. The streamlines snap- 
shots show multi-bubble vortex breakdown evolution, 
convection, merging and shedding. The time-accurate 
integration was carried out up to t = 200 and the so- 
lution showed periodic multi-frequency cycles of vortex- 
breakdown bubbles 21 . An example of the merging of 
vortex breakdown bubbles of same sign of vorticity is 
shown at t = 17. An example of convection and shed- 
ding of vortex breakdown bubbles is shown at t = 25. 
Comparing the streamlines solutions at t = 25 and t = 
89, it is seen that the solutions are almost the same which 
conclusively show that the breakdown process is periodic. 
The Mach-contours show the dynamics of inlet shock mo- 
tion. In the time range of t = 3-8, the inlet shock moves 
upstream toward the inlet and its central portion exists 
outside the inlet section at t = 8. In the time range of 
t = 8*25, the inlet shock moves downstream with corre- 
sponding evolution, convection, merging and shedding of 
breakdown bubbles. In the time range of t = 25-45, the 
inlet shock maintains its motion in the downstream di- 
rection at a slower rate than before, while another shock, 
which is downstream of the inlet shock, appears and also 
moves in the downstream direction. The evolution, con- 
vection and shedding slowly continues until t = 66. In 
the time range of t = 66-78, the downstream shock disap- 
pears and a large vortex-breakdown bubble appears and 
moves upstream. This motion of the bubble is accom- 
panied by upstream motion of the inlet shock (t = 78). 
Later the inlet shock again moves in the downstream di- 
rection and the process is repeated. An animation movie 
has been produced for the total dimensionless time of t 
= 200. Figure 3 shows a blow-up of two snapshots of 
the streamlines solutions at t = 84 and t = 87. At t = 84, 
we recognize five vortex breakdown bubbles which spa- 
tially alternate in their sign of vortex strength. It is also 
noticed that six stagnation points exist at the axis. At t 
= 87, we recognize seven vortex breakdown bubbles and 
seven stagnation points. The figure shows the merging of 
two bubbles of same vorticity sign. 

1.2. Riemann Invariant Boundary Conditions 
with p b = pcc 

In this case, the back pressure at the subsonic points 
of the duct exist, p*, is specified to be equal to poc and 
the other four variables are extrapolated form the inte- 
rior cell centers. The computations have repeated on the 
same grid and for the same flow conditions as that of the 
previous case. Figure 4 shows snapshots of the stream- 
lines and Mach contours of the solution. Comparing the 


present solution with the previous case (Fig. 2), it is seen 
that the two solutions are the same until t = 35. There- 
after, for t > 35, the inlet shock continuously moves in the 
downstream direction with the vortex-breakdown bubbles 
moving ahead of the shock. The shock and vortex bubbles 
are shed and disappears from the duct at advanced levels 
of time. The breakdown mode is termed as “a transient 
multi-bubble vortex breakdown.” The reason behind dis- 
appearance of the shock-vortex-breakdown-bubble sys- 
tem is that the back pressure is so low that it cannot sup- 
port the inlet shock keeping it in the inlet region. More- 
over, the Riemann-invariant type conditions at subsonic 
points let the downstream effects propagate upstream as 
time increases. 

1.3. Riemann Invariant Boundary Conditions 

with p = 2poc 

In this case, the back pressure at the subsonic points 
of the duct exit is specified to be p* = 2poc and the 
other four variables are extrapolated from the interior cell 
centers. Figure 5 shows snapshots of the streamlines and 
Mach contours of the solution. Comparing the present 
solution with the solution of the first case (Fig. 2), it is 
seen that the two solutions are similar with the exception 
that the present solution lags that of the first case in 
time. The reason behind this behavior is that the back 
pressure p* is larger than that of the first case. Moreover, 
the Riemann invariant type conditions at subsonic points 
let the downstream effects propagate upstream as time 
increases. The existence of the large back pressure which 
is felt upstream supports the inlet shock and keeps it in 
the inlet region. 

1.4. Riemann Invariant Boundary Conditions 

with |f = constant 

In this case, the back pressure is obtained from the 
condition |f = constant at the subsonic points of the duct 
exit The other flow variables are extrapolated from the 
interior cell centers. Figure 6 shows snapshots of the 
streamlines and Mach contours of the solution. Compar- 
ing the present solution with the solution of the first case 
(Fig. 2), it is seen that the two solutions are similar until 
t = 22. Thereafter, for t > 22, the inlet shock continu- 
ously moves in the downstream direction with the vortex- 
breakdown bubbles moving ahead of the shock. Again, 
as in the case of Fig. 4, the shock and vortex bubbles 
are shed and disappear from the duct at advanced levels 
of time. The breakdown is termed as “a transient multi- 
bubble vortex breakdown.” The reason behind disappear- 
ance of the shock-vortex-breakdown-bubble system is that 
the back pressure obtained from |f = constant condition 
is so low that it cannot support the inlet shock keeping it 
in the inlet region. Moreover, the Riemann-invariant type 
conditions at subsonic points let the downstream effects 
propagate upstream as time increases. 


1.5. Placing a Disk at the Exit with r = 0.333 

In (his case, a circular solid disk of radius r * 0.333 is 
placed at the duct exit and solid-boundary conditions are 
applied on the disk surface. For the remaining portion 
of the exit, the boundary conditions are obtained by 
extrapolation from the interior cell centers. Figure 7 
shows snapshots of the streamlines and Mach contours 
of the solution. It is noticed from the Mach contours 
that most of the exit points are subsonic. The streamlines 
show that the shock-vortex-breakown-bubble system first 
appears behind the duct exit. Thereafter, they move 
upstreams until the inlet shock moves outside of the inlet. 
The vortex -breakdown bubbles are then locked between 
the shock outside of the inlet and the circular disk at the 
exit 

II. Supersonic Swirling Jet from a Nozzle 

Figure 8 shows the computational domain and a typ- 
ical grid for this external flow case. The dimensions of 
the computational domain is 7 x3.5 in the axial and r adial 
directions, respectively. The grid is 221x51x2 points in 
the axial, radial and tangential directions, respectively. 
The grid is clustered at the nozzle exit (x = 0.0) and at 
the domain centerline. The dimensionless nozzle radius 
is unity, where a supersonic jet of M x = 3.0 is issued, 
and outside the nozzle another jet is issued at M x = 2.0. 
The freestream Reynolds number is 296,000. 

Figure 9 shows the inlet flow profiles of the ax- 
ial velocity, swirl velocity, radial velocity, pressure and 
density, which are taken from the experimental data 
of Ref. 12. The initial profiles are used as quasiaxi- 
symmetric profiles for the present computations. On the 
cylindrical boundary (side wall) of the flow at r = 3.5, 
freestream conditions are imposed corresponding to M x 
= 2.0. The initial conditions in the computational domain 
are also taken as those corresponding to the freestream 
conditions at M x = 2.0. The problem is solved using 
two types of exit boundary conditions at x = 7.0; first ex- 
trapolation of all five variables from the interior cell cen- 
ter and second using the Riemnn-invariant-type boundary 
conditions. 

II.l. Extrapolation from Interior Cell Centers 

Figure 10 shows snapshots of streamlines and Mach 
contours of the solution. The streamlines show multi- 
bubble breakdown at the early levels of time. These bub- 
bles develop due to the shock system formed at the nozzle 
exit in the vicinity of the centerline. It is noticed that a 
strong portion of the shock exists at the centerline which 
splits into two oblique shocks, one is a weak shock and 
the other is a strong shock. Behind the strong shock, the 
vortex breakdown bubbles exist Thereafter, for t > 5, the 
oblique shocks move slowly in the downstream direction 
and breakdown the vortex bubble stays in its place At t 
> 55, the shock system moves upstream and so does the 
vortex breakdown bubble. The slow motion of the shock 


system and the vortex breakdown bubble continues back 
and forth between these two locations. No vortex shed- 
ding has been captured during the computations of this 
case. It is also noticed that most of the exit points are 
continuously supersonic and hence no downstream effects 
exist with the exception of a very thin-layer around the 
centerline. 

II.2. Riemann Invariant Boundary Conditions 

Next, the boundary conditions at the exit are replaced 
by using the Riemann-invariant-type boundary conditions 
with p t = p x at the subsonic points. Figure 1 1 shows 
snapshots of the streamlines and Mach contours of the 
solution. By comparing the present solution with the 
previous case of Fig. 10, we see that there is very slight 
effect of the present boundary condition on the solution. 
This is understood since the subsonic region at the exit is 
very small and moreover, the exit boundary is far from 
the nozzle exit 

Figure 12 shows a blow-up of the Mach contours at t 
= 55 for the flow case of Fig. 10. The shock system near 
the nozzle exit is clearly seen. 

Concluding Remarks 

The unsteady, compressible full Navier-Stokes equa- 
tions are used to study the critical effects of the down- 
stream boundary conditions on the supersonic vortex 
breakdown. In the present study, two supersonic swirling 
flow cases are considered. The first one is for a super- 
sonic swirling flow in a configured circular duct, where 
four types of exit boundary conditions are used. The sec- 
ond one is for a supersonic swirling jet issued from a 
nozzle into another supersonic jet of lower Mach num- 
ber. This flow case is a computational simulation to the 
experimental case of Ref. 12. 

In the first flow case, we have shown that the 
downstream exit boundary conditions produce substantial 
changes in the vortex breakdown modes, particularly at 
advanced time levels. This is due to the effect of the sub- 
sonic downstream points on the flow upstream. In the sec- 
ond flow case, we have shown that the downstream exit 
boundary conditions do not produce substantial changes 
in the vortex breakdown modes. This is due to the very 
thin subsonic region at the exit and the large distance of 
the downstream exit from the nozzle exit 
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Figure 3. Blow-up of streamlines of periodic multibubble, multi- 
frequency vortex breakdown, extrapolation from interior 


B.C. 




















Figure 5. Streamlines and Mach contours for a swirling 
flow with periodic multibubble, multifre- 
quency vortex breakdown, P t = 2 Poo, Rie- 
mann invariant B.C. 






















Figure 7. Streamlines and Mach contours for a swirling flow 
with quasi-steady multibubble vortex breakdown, down- 


stream disk of r = 0.333. 
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re 10. Streamlines and Mach contours for super- 
sonic swirling jet from a nozzle wuh low- 
frequency almost single bubble vortex break 





Figure 11. Streamlines and Mach contours for super- 
sonic swiriing jet from a nozzle with low 
frequency almost single bubble vortex break- 
down, Riemann invariant B.C. 
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Figure 12. Blow-up of Mach contours at t = 55 for supersonic 
swirling jet from a nozzle with low frequency almost 
single bubble vortex breakdown, Riemann invariant 
B.C. 





